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ABSTRACT 

Exoplanetary systems are found not only among single stars, but also binaries of 
widely varying parameters. Binaries with separations of 100-1000 au are prevalent in 
the Solar neighborhood; at these separations planet formation around a binary member 
may largely proceed as if around a single star. During the early dynamical evolution of 
a planetary system, planet-planet scattering can eject planets from a star's grasp. In a 
binary, the motion of a planet ejected from one star has effectively entered a restricted 
three-body system consisting of itself and the two stars, and the equations of motion 
of the three body problem will apply as long as the ejected planet remains far from 
the remaining planets. Depending on its energy, escape from the binary as a whole 
may be impossible or delayed until the three-body approximation breaks down, and 
further close interactions with its planetary siblings boost its energy when it passes 
close to its parent star. Until then this planet may be able to transition from the 
space around one star to the other, and chaotically 'bounce' back and forth. In this 
paper we directly simulate scattering planetary systems that are around one member 
of a circular binary, and quantify the frequency of bouncing in scattered planets. We 
find that a great majority (70 to 85 per cent) of ejected planets will pass at least 
once through the space of it's host's binary companion, and depending on the binary 
parameters about 45 to 75 per cent will begin bouncing. The time spent bouncing is 
roughly log-normally distributed with a peak at about 10 4 years, with only a small 
percentage bouncing for more than a Myr. This process may perturb and possibly 
incite instability among existing planets around the companion star. In rare cases, the 
presence of multiple planets orbiting both stars may cause post-bouncing capture or 
planetary swapping. 

Key words: binaries:general-planet-star interactions-planets and satellites: dynam- 
ical evolution and stability-scattering 



1 INTRODUCTION 

Binary star systems which contain exoplanets vary dramat- 
ically. From the fir st discoveries of ex oplanets orbiting stars 
in binary systems (|Butler et al.lll997l) to a pul sar planet with 
a stellar companion (|Sigurdsson et all [2003) to the recent 
circumbinary "Tat ooine"-like planet s announced by the Ke- 
pler mission team (|Dovle et alJbOlllJWelsh et al.ll2012r ). ex- 
oplanets have been observed to exist in an intriguing variety 
of multiple-star systems. A few years after the discovery of 
the fi rst multi-planet sys t em around a main sequ ence star, v 
And l|Butler et aiJll999F ). lLowrance et"al] |2002T ) discovered 
a stellar companio n to v And at ~ 750 AU. Subsequently, 
ICuriel et alj (|201ll) announced a fourth planet orbiting the 
primary (see also McArthur et al.l 2010). The orbital archi- 
tecture of this system may be explained by secular interac- 
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tions sub sequent to planet -planet scattering and planetary 
ejection (|Ford et al.ll2005[). and /or pe rturbations from the 
stellar companion (|Barnes et al.ll201lh . 

Generally, systems with multiple stellar and planetary 
components present rich dynamical laboratories which could 
inform planetary formation and evolution theories. Indeed, 
iDesidera fc Barbieril (|20Q7m estimate that binary compan- 
ions within just 100 — 300 AU are likely to play a role 
in the formation and evolution of planetary companions. 
In the Solar ne i ghbou rhood, such binaries are prevalent; 
Rash avan et ail (|2010h . in the most complete survey our lo- 
cal volume within 25 pc, find approximately half of solar- 
type stars in multiple systems, with a broad peak in the 
(roughly lognormal) separation distribution stretching from 
a few to a few thousand au. At separations above a few tens 
of au, these local multiple systems have the same planet 
hosting frequency as single stars, setting the stage for po- 
tential interac tions between s t ellar and planetary dynamics. 
An example is lMarzari et al.l (|2005l h who studied planetary 



© XXXX RAS 



2 Moeckel & Veras 



systems undergoing planet-planet scattering from around a 
star with a binary companion. While the focus of that work 
was on the orbits of the unejected planets, they found that 
the binary companion participated directly in the planet 
scattering, and a few per cent of the ejectees collided with 
the binary partner; interaction between the binary star and 
the planets of the most direct sort. 

Here, we investigate planetary transfer between stellar 
compan ions subsequent to p lanet-planet scattering. In con- 
trast to M arzari et al l (|2005l ). we are most interested in the 
planets that are ejected from their host star rather than the 
planetary system left behind. Our motivation is best illus- 
trated by considering the motion of test particles in the pla- 
nar circular restricted three body problem (planar CR3BP). 
The planar CR3BP, whose early form was establishe d by 
Euler and Lagrange in 1772 (seeH 

arrow- Greenlll997l ). de- 
scribes the motion of a massless particle subjected to the 
force of two massive bodies moving in circular orbits such 
that all bodies are coplanar.The CR3BP (in its planar form 
or its extension into thre e dimensions, th e 3D CR3BP, e.g. 
ISzebehelv fc Peters! Il967l ; iMarchall Il990h is well suited to 
describing the motion of an extrasolar planet orbiting one 
or both stars of a circular binar y system; even the most 
massive planets, at about 13 Mj (|Burrows et al.|[T997l ). are 
< 10 per cent as massive as the lowest-known exoplanet host 
star mass (as of 31 January, 2012, 0.13 M Q for KOI-69lfl 
However, in the vast majority of cases, the planet/star mass 
ratio is much less than 1 per cent. Hence, in binary systems 
with multiple planets, the motion of a planet will be dic- 
tated by the CR3BP as long as the perturbations from the 
other planets are negligible. This situation arises when one 
planet is sufficiently far away from the others, which may 
occur subsequent to planet-planet scattering. 

The position and velocity of this planet, as well as the 
masses and separations of the two stars, d etermin e what 
regimes of motion are possible for the planet, [jacobil (|l836l ) 
showed that these quantities are related by a constant of the 
motion, C, such that 

C 2 + !/ 2 + 2f^ + ^^ 2 -y 2 -i 2 , (1) 
\ r v r s J 

where the position (x, y, z) and velocity (x, y, z) of the planet 
are computed with respect to the origin of a frame rotating 
with a mean motion of unity about the center of mass of 
the stars such that /i p = GM P , fi s = GM S , y, p -\- fj, s — 1 
and where r v and r s represent the distance of t he pl a net to 
the primary and secondary, respectively. Later, iHilll (|l8T8l ) 
linked given values of C with regions of space which are 
allowed or forbidden. 

Figure [1] shows examples of these regions in the planar 
case, i.e. z = z — in equation [T] Suppose a massless planet 
is initially orbiting the primary. Then, the trajectory of the 
planet is given by the continuous line, the primary and the 
secondary are represented by the large red and blue dots 
on the left and right respectively, and the gray areas show 
regions that the planet cannot access. The relative masses 
of the stars here are /^s/m p = 0.3. In the C — 5 plot, the 
planet is restricted to orbiting the primary. In the C = 3.7 
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Figure 1. Curves of zero velocity in the dimensionless synodic 
frame of a circular binary system. In this frame the stars (red and 
blue dots) are fixed in space. The gray regions show inaccessible 
regions for test particles with the labeled value of the Jacobian in- 
tegral. The continuous lines show example orbits at these energies, 
each originating in the space around the more massive star. At 
even lower values of the Jacobian constant, the forbidden regions 
contract around the L4 and L5 Lagrange points before vanishing. 

plot, the planet may orbit either the primary or secondary, 
and further, freely travel in between each star's sphere of 
influence via the connecting region centered on the stars' 
Li Lagrange point. The bottom two plots illustrate different 
ways in which the planet can escape the system; the C = 3.1 
plot demonstrates how the planet can orbit both stars before 
escaping. While this illustration is in the framework of the 
planar CR3 BP, such transiting orbits exist in the 3D case 
as well (e.g. lGomez et alJliool i. 

A planet forming around one star of a widelj0 sepa- 
rated binary will be in the limit shown in the upper left 
panel of Figure [1] although the C value of this planet will 
constantly change if there are other planets around one or 
both stars. However, this change is significant only during 
close approaches with the other planets. Planet-planet scat- 
tering may cause a sudden decrease in C, enough to open up 
the gap around L\ seen in the C = 3.7 plot but not enough 
to enable escape the system through the exterior L2 and L3 
Lagrange points, as seen in the C = 3.5 and C = 3.1 plots. 
In this case, the planet may 'bounce' between the stars. If 
the planet is scattered far enough away from the other plan- 
ets, then C might maintain a value that allows this behavior 
to continue on appreciable timescales. Here, we explore this 
possibility, and characterize the prospects and consequences 
for planetary escape and capture that is internal to a binary 
stellar system. 

Because planets are not massless, and we wish to treat 
planet-planet scattering accurately, we cannot rely solely on 

2 'Widely' for these purposes is loosely defined as a system with 
semi-major axes such that denary 3> a-pianets, so that the planets 
are initially orbiting a single star. 
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e=0.999 binary Pythagorean problem Kozai oscillations 




orbits time / Myr 



Figure 2. Tests of the integrator used in this paper. Left: Integration of a planet with eccentricity e = 0.999 for 1000 orbital periods. 
The errors in energy and angular momentum are small and periodic, with no systematic drift. There is, however, a small linear drift in 
the longitude of perihelion. Middle: Paths of the stars in the Pythagorean problem t hrough t = 80. In each panel the positions of the 
bodies at the end of the time interval are shown with dots. Right: A test presented in lNaoz et al a system with an inner binary 

consisting of a 1 Mq star and a 1 Mj planet on an a = 6 au, e = 0.001, i = 64.7° orbit, with a 40 Mj companion on an a = 100 
au, e = 0.6, i = 0.3° orbit. The longitude of pericenter of both companions are initially zero. The inclination and eccentricity of the 
inner planet ar e plotted as it und ergoes Kozai oscillations, including switches from prograde to retrograde orbits. The agreement to the 
integrations in lNaoz et al.l (1201 lh is excellent. 



the analytical structure of the planar or 3D CR3BP for our 
explorations. In particular, planets often reach ejection ve- 
locities via a series of encounters that send them progres- 
sively farther from their host star. This rules out introduc- 
ing test particles with a spectrum of energies in the CR3BP. 
If ejections were dominantly the result of single encounters, 
the CR3BP would be more readily and directly applicable, 
although even then we would need to introduce a stochastic 
forcing effect within a certain distance of the host star to ac- 
count for possible interactions with the remaining planets. 
We therefore turn to direct numerical simulation. Because 
A^-body codes designed to study planetary systems often rely 
on the presence of one massive central object that dominates 
the motion of the other bodies, we constructed a simple and 
flexible integrator to investigate planetary bouncing. This 
code and its tests are described in Section[2] In Section[3] we 
describe ensembles of simulations with the code, and show 
that bouncing is the dominant behaviour of planets leaving 
a binary system via planet-planet scattering. We discuss the 
results and future extensions of this work in Section 3] and 
then briefly conclude in Section 



2 THE CODE 

The ./V-body code we use is built upon a basic fourth-order 
Hermite integrator, written by P. Hu10. The straightforward 
structure of the code makes it easy to adapt for different 
purposes, although this comes at the sacrifice of fine-tuning 
that could speed up the code if it were tailored to specific 
problems. Our modifications and tests are detailed here. 



2.1 Integration Algorthim 

In its implicit form, the Hermite method is time symmetric 
and exhibits no drift in energy when integrating periodic or- 
bits. In practice, an explicit approximation to the implicit 
method is usually used, following a PEC (predict-evaluate- 
correct) cycle. This approach is commonly applied to stel- 
lar dynamics probl ems, where supplementary met hods (e.g. 
KS regularisation, iKustaanheimo fc St icfcl 1965) are used 
to avoid systematic errors in long lived systems with pe- 
riodic motion. In order for the method to be suitable for 
studying planetary systems, the implicit version of the Her- 
mite integrator must be used. To achieve this we make two 
modifications to the base code. First, we iterate over the 
evaluation of the force and its time derivative and the correc- 
tion to the predicted values, converting to a P(EC) n method 
l|Kokubo et al.f l998). This iteration process converges to the 



3 That code is available at http://www.artcompsci.org 
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implicit Hermite solution, at the expense of multiple force 
evaluations. In practice, we use a single iteration, n = 2. 

With fixed timesteps, the P(EC) n method would yield 
excellent energy conservation. When timesteps are allowed 
to be variable, the timesteps must be chosen to be symmet- 
rical functions of the system's state at the beginning and 
end of the timestep in order to retain the energy conserving 
features of the Hermite scheme. When all the particles share 
the same timestep, as in our code, t his is straightforw ard to 
implement, following the method of lHut et al.l ((l995). This 
approach requires iterations over the timestep choice, and we 
again use one iteration here. A single timestep then involves 
one iteration for the P(EC) n integrator, wrapped inside an- 
other iteration for the timestep determination. The result is 
a robust and flexible integrator for small- iV systems with no 
preferred dominant force or geometry. 

The timestep requested by each pair of particles i,j 
with relative separation ry , velocity , and acceleration aij 
is proportional to the minimum of their estimated free fall 
and collision times, the latter under unaccelerated motion: 
Stij = T)mm[(aij /rij) 1 / 2 ,rij /vij]. The system timestep for 
a given particle state is then the minimum of this quantity 
over all particle pairs; this estimate is used in the timestep 
symmetry iteration. In the simulations presented here we 
set the control parameter rj = 0.02, yielding typical rela- 
tive energy errors of the order 10 -8 , and relative angular 
momentum errors of the order 10" 11 . 



2.2 Tests 

In order to verify the integrator's capabilities to follow plan- 
etary dynamics, we have performed several tests. We show 
the results in Figure [5] In the first, we integrate a binary 
system with eccentricity e = 0.999 for 1000 orbits. We plot 
the semi-major axis, eccentricity, and longitude of perihelion 
over the run of the integration. The numerical errors are pe- 
riodic with minimal systematic drift, with the exception of 
a small linear error in the longitude of perihelion. Similarly 
good results are obtained with less extreme eccentricities. 
The example shown has a secondary to primary mass ratio 
of 0.1, and similar results are obtained for ratios from unity 

to approximately zero. 

The second test is of the Pythagorean problem jBurraul 

1913). With G = 1, three bodies of mass 3, 4, and 5 are ini- 
tially at rest at the corners of a Pythagorean right triangle, 
each body at the corner opposite the side of the triangle cor- 
responding to its mass. For a detailed ac count of the complex 
three body dance that ensues, refer to Szcbchcl v fe Peters! 
(1 19671 ). the first authors to resolve the system to its conclu- 
sion: after one body is ejected from the system, while the 
other two recoil as a binary. Those authors used a regulari- 
sation scheme to deal with the difficulties of the many close 
encounters in the system's evolution. Using the parameter 
for the timestep control that we use throughout this paper, 
our code integrates the system through virtually the same 
trajectory as Szebeheley & Peters find, with relative energy 
error of the orde r 10 -9 . Following Szebehely & Peters and 
iHut et ail (|l995l ). we integrate the system with velocities 
reversed from t = 62 back to t = 0. Achieving similar re- 
sults as those authors, we recover the initial conditions with 
errors entering at the third decimal place. 

Finally, we integrate a system undergoing Kozai oscil- 
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Figure 3. Semi-major axes and eccentricities of planets in a 250 
au, equal-mass binary. The orbital elements are plotted relative 
to whichever star's space the planet is found in at any time, as 
defined in the text. Top panels are orbits around the host, bottom 
panels are around the companion. 



lations (lKozai|[l962f ). leading to flipping of the inclination 
from prograde to retrograde. The setup is as follows: a 1 
Mq primary is orbited by a 1 Mj planet with a = 6 au, 
e = 0.001, and i — 64.7°. There is also a 40 Mj companion 
with a — 100 au, e = 0.6, inclined 65° with respect to the 
inner planet. As the eccentricity and inclination of the com- 
panions oscillate, the inner planet oscillates between pro- 
grade and retrogra de orbits. This is the same system pre- 
sented in figure 3 of lNaoz et al.l (|201lf ). to which our results 
can be compared. Those authors integrated the system us- 
ing a Burlisch-Stoer integrator, and the agreement with that 
result is excellent. 



3 THE SIMULATIONS 

We performed several sets of 300 simulations each, exploring 
a limited but illustrative range of binary parameters. Each 
simulation is run for 10 Myr. 



3.1 Initial Conditions 

We follow a stan dard planetary setup for scatterin g sim- 
ulations, following iMarzari fe Weidenschillind (|2002f ). Each 
adjacent pair of planets is separated by K mutual Hill radii, 
so that the planet i and its outward neighbour j have semi 
major axes Oj — at + KRmj, with 

, 1/3 



rxii + rrij 



di + a,- 



and the chosen val ue of K determines th e insta- 
bility timesca le ( Marzari fe Weidenschillingl 120021 : 

IChatteriee et all 12008 ). Planetary masses are drawn 
from a power law mass function with f(m) oc m~ , 
between 0.3 and 5.0 Mj. We take K = 4 and an innermost 
semi-major axis of 3 au, which together with our mass 
spectrum yields an instability timescale of the order 10 4 
yr, although the possibility of long-term stable systems 
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Table 1. Parameters and scattering results for the simulated binary systems. All binaries have zero eccentricity. The majority of scattering 
planetary systems experience a bouncing planet. 



a j a,u 


Mt. /Mn 




•A/" scatter 


1 ' switch 




ejections 






collisions 














Ntotal 


^switch 


^bounce 


Ntotal 


"^switch 


"^bounce 


250 
250 
250 


1.0 
1.0 

0.3 


1.0 
0.3 
1.0 


249 
231 
209 


221 
199 
200 


236 
220 
226 


204 (86±7%) 
192 (S7±l%) 
185 (82±6%) 


182 (77±|%) 
138 (63±|%) 
174 (77±6%) 


84 
67 
43 


14 (17+^%) 
4 C6.0±*;S%) 
10 (23±^°%) 


(0l^ 2 %) 
(0 +2 7 %) 
3 (7.0l|:|%) 


1000 
1000 
1000 


1.0 
1.0 

0.3 


1.0 
0.3 
1.0 


235 
227 
165 


185 
157 
116 


215 
207 
162 


180 (Utl%) 
156 (75±l%) 
112 (69±7%) 


133 (62^%) 
73 (35±|%) 
88 (54±6%) 


69 
73 
21 


1 (1.4±i;*%) 

(o+g- 5 %) 

2 (8.3j*l%) 


(0± 2 2 %) 
(0 +2 5 %) 

o 



Mj, and M c are the masses of the host and the companion stars. Mscatter is the number of systems that scattered, out of 300 sets of initial 
conditions. N gw i tch is the number of planets that cross the L\ threshold at least once. In the 'ejections' and 'collisions' columns, n sw i tch 
and ni, ourlce are the number that cross L\ at least once and more than once, respectively; the percentages include Poisson confidence 
limits corresponding to a la Gaussian interval. 



250 au binary 




time around host / yr time around host / yr time around host / yr 

1000 au binary 




time around host / yr time around host / yr time around host / yr 

Figure 4. Scatter plots of the time spent in the space around each star for all planets that bounced across the boundary between them. 
The zero point for the clock is the moment the planet first crosses over the threshold; all the time spent around host prior to being 
scattered across to the companion is not counted. The slight gridded appearance of points at small times is due to our time resolution 
for these plots, which is 250 years. The diagonal gray line shows equal time spent around each star. The error bars at the bottom and 
left show the median of the distribution, and the 16th and 84th centiles. 
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exists with this setup. Eccentricities are initially zero, 
inclinations are uniformly distributed between 0° and 1°, 
and all the phase angles are randomized. Changing the 
planets' eccentricity and mutual inclinations will alter the 
inst ability timescale, but ot herwise yield similar outcomes 
(e.g. lJuric fc Tremaindl2008h . 

In this initial work, the binary system initially has zero 
eccentricity, and has semi-major axes of 250 or 1000 au. 
Small changes to these values are possible as the system 
evolves and planets are ejected or collide with the stars. 
The semi-major axis values are chosen to lie within the peak 
of th e local binary separation distribution (IRa ghavan et al.l 
2010), but large enough that the planet formation process 
around a single star is probably not drastically affected by 
the presence of its binary companion. The masses of the 
planet hosting star and its binary partner (henceforth the 
'host' and 'companion', respectively) are taken to be (1.0, 
1.0), (1.0, 0.3), and (0.3, 1.0) Mq. We assume coplanarity 
between the the orbital plane of the binary and the plan- 
etary system; this should have only a secondary effect on 
the results, since planetary ejections from scattering do not 
occur strictly in the plane of the planetary system. 

A planet is considered to be ejected if its distance from 
the binary center of mass exceeds 5 x 10 4 au, a factor of 
~ 2 less than a typical distance at which galactic tides can 
cause escape over a typical main sequence lifetime (|Tremainei 
1 19931 ) but well beyond the region in which an appreciable 
number of planets may be scattered onto stable wide orbits 
I Veras et al.ll2009r i. Collisions between bodies occur when 
their radii are detected to overlap; the planets are all 1.3 
Jupiter radii, and the stars are 5 Solar radii. Using this 
large stellar radius enables us to avoid the regime where 
tidal interactions overly affect the dynamics. Mergers be- 
tween bodies are momentum and mass conserving. 

3.2 Analysis and Characterisation 

To assign a planet at any instant to one of the two stars, 
we use the following simple scheme: first, rotate the system 
about the binary's center of mass so that the stars are on 
the x-axis, and the binary's angular momentum is in the z 
direction. The location of the binary's L\ Lagrange point is 
determined, and this is compared to the positions along the 
i-axis of the planets; whichever side of the L\ divide a planet 
lies on is its instantaneous host. Planets with distances from 
the binary center of mass greater than twice the separation 
of the binary components are tagged so that a planet on a 
long looping excursion, or one that has been ejected but not 
yet removed from the calculation, is not counted as being in 
orbit about either star. 

While quite simple, this scheme is motivated by the 
symmetry of the allowed regions of space in the CR3BP 
(see Figure [1]). After subtracting off the orbital motion of 
the binary and appropriately nondimensionalising all veloc- 
ities, masses, and distances, we also track an approximation 
to the Jacobian integral of each planet in the 3D CR3BP. 
In doing this we treat the host and its associated planets 
as a single body at their barycenter. While this serves to 
estimate when a planet has sufficient energy to transition 
between topologies of allowed orbits, it is not exact: the 
non-zero mass of the planets and small eccentricities in the 
binary orbit induced when a planet is ejected restrict this to 



an approximation. With the history of each planet's stellar 
residence recorded, we tally the number of times each planet 
crosses the L\ threshold, as well as the amount of time spent 
around each star. When a planet crosses the Li threshold 
a single time, we say that it has 'switched'; if it crosses the 
threshold repeatedly, it has begun to 'bounce'. 

The main questions we wish to answer here are: what 
fraction of scattered planets transition to the space sur- 
rounding the companion star, and how much time do they 
spend around the companion, or bouncing between the 
stars? 

3.3 250 au binary results 

We turn first to the results for a 250 au binary. Table 
12.21 provides a summary of the outcomes from our sim- 
ulations. For the systems with 1 Mq hosts the number 
of systems that go unstable and scatter (Abutter) dur- 
ing the 10 Myr of our simulations is about 230-250, i.e. 
15 to 20 per cent of the runs have not yet scatterecjf]. 
The systems with a 0.3 Mq host seem to be more stable, 
with about 30 per cent unscattered. Some difference is ex- 
pected; t he instability timescal e is dependent on the mass 
funct ion ( Chamber s et al. 1 1 19961 : IM arzari fc Weidensch illing 
2002). Keeping the planetary mass function fixed while re- 
ducing the host mass is equivalent to shifting the planets to 
higher masses. The sense of the change in instability times 
we see is the same as in previous work: relatively more mas- 
sive planets are stable for longer times when spaced as fixed 
multiples of their mutual Hill radi{f|. 

Looking at the number of planets that switch, we see 
that a large fraction of scattered planets will pass through 
the companion star's space; the number of switching planets 
N SW itch is about 200 for all three binary setups. The over- 
whelming contribution to this category comes from plan- 
ets that end up ejected from the binary and eventually re- 
moved from the calculation, shown the 'ejec- 
tions' columns. For instance, 204 of the 221 planets that 
switch in the equal-mass binary runs end up ejected, with 
similar fractions for the other two mass ratios. Crossing the 
L\ threshold is symptomatic of a planet that is leaving the 
system. We illustrate one such system in Figure [3] which 
plots the semi-major axes and eccentricities of the planets 
with respect to whichever star's airspace they are in at the 
time. The scattering in this case takes place very early, and 
the planet shown by the blue line begins to transition back 
and forth. After an uninterrupted period of about 0.5 Myr 
around the companion, the final few bounces back to the 
host bring it into further interaction with the other planets, 
and its next trip through the companion's space ends with 
escape from the system. 

Looking further at the ejected planets, note that the 
majority (about 85 per cent) cross over to the companion 

4 This fraction is consist ent with the results of an id entical plan- 
etary setup presented in iMoeckel fc Armitagei i2012h , integrated 
wi th MERCURY dChamberslll999h . 

5 IChambers et aU l[l996) suggest that the instability timescale in 
a three planet system should vary independently of the planets' 
mass function when the planets' spacing is scaled as the one- 
fourth power of their masses, rather than the one-third power in 
the Hill radius spacing scheme used here. 
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on their way out of the system: ejection from the host star 
occurs predominantly through the space around L\ rather 
than the exterior Lagrange point. Furthermore, a majority 
of ejected planets (about 65 to 75 per cent) not only switch 
but bounce, traversing between stars multiple times. Esca- 
pers are not simply passing through the companion star's 
space, but rather orbiting there for some time before cross- 
ing back to the host, and so onward chaotically until they 
leave the system or interact further with the rest of the plan- 
ets. Planets that ultimately collide with something are less 
interesting from a switching standpoint; very few collision 
partners are involved in intrabinary excursions. 

The answer to the question of what percentage of 
ejected planets cross over to orbit the companion is evi- 
dently an interestingly large number: something like 85 per 
cent of escapers will switch. In Figure [4] we address the time 
they spend in the space around each star by plotting the the 
time around the host versus the time around the compan- 
ion for all stars that make the jump and return, with the 
clock starting from the moment they first pass through the 
L\ connection-i.e. we do not count time spent around the 
host prior to scattering. The diagonal line in the plots shows 
equal time around each star; the points cluster around this 
line in a clear correlation. The point lying near 1 Myr around 
each star in the equal mass case is the system shown in Fig- 
ure [3] A more typical bouncing system is shorter lived, with 
fewer long term residencies about either star. There is a bias 
towards spending time around the more massive star in the 
unequal mass binary cases, but this is not a strong effect. 
There are a few outliers in each case with large times around 
the host. These are all cases where a planet made excur- 
sions to orbit around the companion, and at some point fur- 
ther interactions with the planets around the host changed 
the Jacobian integral of the wandering planet, re-isolating 
it around the host. These adjustments to the Jacobian inte- 
gral are not possible when the planet is in orbit around the 
companion. 
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Figure 5. The distribution of times that bouncing planets spend 
around each star in years (top panel) and binary orbital periods 
(bottom panel). Time around the host star are shown with solid 
lines, and time around the companion as dashed lines. The medi- 
ans of the distributions are shown with arrows at the top of each 
plot. 
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3.4 1000 au binary results 

Overall, the results for the 1000 au binary are quite similar 
to the 250 au case. This is largely expected; the spacing of 
the planets, set at a fixed inner radius of 3 au, together with 
the physical radii of the bodies are the physical quantities 
that generate the ejection velocities of the planets. The ratio 
of the planets' semi-major axes around the host to the binary 
separation is very small for both of the binary sizes, and the 
orbits onto which the planets are launched when they begin 
bouncing should originate from similar points in space when 
the binary separation is taken as the length scale. The only 
remaining explanation for changes in the outcomes is the 
relation between the ejection velocities, a physical quantity 
set by the planetary dynamics, and the Jacobian integral. A 
planet with a given physical velocity and position relative to 
its host will have a lower Jacobian integral in a more widely 
spaced binary, since the dimensionless velocity is tied to the 
binary period. 

What differences there are in the outcomes compared 
to the 250 au binaries are consistent with this picture. The 
percentage of ejected planets that switch are essentially the 
same for both binary separations, but the bouncing per- 
centage is lower for the 1000 au binaries. A planet that at- 




1000 au binary 

250 au binary 



2.5 3.0 3.5 

Jacobian integral at switch 



4.0 



Figure 6. The bottom panel shows the distribution of the Ja- 
cobian integral C at the moment a planet first switches to the 
companion for the equal mass binaries. The top panel shows, in 
grey, the restricted regions of space in the synodic coordinate sys- 
tem of the CR3BP for three values: the median values of the 1000 
and 250 au distributions, and the 75 centile of the 250 au runs. 



tains the same ejection velocity from around its host in the 
1000 au binary will have a more 'wide open' escape through 
the companion's exterior Lagrange point; this easier escape 
plausibly reduces the tendency to bounce back to the host. 

The times spent around each star after bouncing begins 
are longer for the 1000 au systems. If the systems were simple 
rescalings of an identical dimensionless setup, this difference 
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should be eliminated by rescaling the time to the binary or- 
bital period. We test this with the equal mass binary cases. 
In Figure [5] we show the distributions of time spent around 
each star, which are roughly log-normal, for both binary sep- 
arations. With time measured in the binary orbital period 
(the lower panel) the histograms become more similar, al- 
though the relative positions of the histogram peaks swap 
places. Once a planet has been launched by the internal dy- 
namics of the planets around the host into the three-body 
system consisting of the binary and the ejected planet, the 
planets bouncing in the 250 au binary are longer lived rela- 
tive to the binary's dynamical timescale. This is again con- 
sistent with the idea that the tighter binary's planets find 
themselves with a higher Jacobian integral, and less readily 
able to escape through the region surrounding the exterior 
Lagrange points. 

We can directly test this idea by examining the distribu- 
tion of the Jacobian integral C at the moment of a planet's 
first switch. In Figure [6] we plot this distribution for the 
equal mass cases for both 1000 au and 250 au binaries. We 
also show the intersection with the binary's orbital plane of 
the zero- velocity surface associated with three values of the 
Jacobian integral: the median value for each binary separa- 
tion, and the 75th centile of the 250 au binary. This latter 
value is about the 95th centile of the 1000 au runs. The crit- 
ical value of the Jacobian integral when the L± connection 
opens is Ci = 4.0; at 62,3 = 3.46 L2 and L3 open up. As 
noted earlier, the calculated value of C is not exact, due to 
small errors introduced by the inconsistencies with the 3D 
CR3BP in our actual setups. Despite these minor limita- 
tions, the difference in the distribution of C is clear; more 
planets around the 250 au binary have 62,3 < C < C%, 
where access to the outside universe through the exterior 
Lagrange points is impossible without further perturbations 
from the remaining planets. In contrast, the great majority 
of switching planets in the 1000 au binary are free to escape 
the system. 

Since we observe a dimensional universe, the longer 
physical timescale associated with the 1000 au binaries is 
of potential interest. In particular there are a handful of 
cases that spend long periods around the companion, either 
switching back and forth repeatedly or in rare cases ending 
up in a quasi-stable orbit about the companion. An exam- 
ple of the latter outcome is shown in Figure [7] In this case a 
planet bounces between stars for several Myr, before settling 
around the companion at about 3 Myr and staying there for 
the remainder of the simulation. We stress that this outcome 
is rare, and is not reliably stable. The semi-major axis of the 
planet around the companion is about 200 au, and the con- 
tinued influence of the host star is evident in the eccentricity, 
which oscillates between about 0.5 and 0.95. 

3.5 Response of a planet around the companion 

The orbits of stars passing through the companion's space, 
even those that are only bouncing for a very short amount 
of time, are characterised by high eccentricities. Of potential 
interest is the effect this may have on planets or discs around 
the companion. We examine this very briefly by conducting 
simulations identical to the 250 au, equal mass binary runs, 
but with the addition of a single planet around the com- 
panion. This planet has a mass of 1 Mj, a semi- major axis 
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Figure 7. Semi-major axes and eccentricities of planets in a 1000 
au, equal-mass binary. The orbital elements are plotted relative 
to whichever star's space the planet is found in at any time, as 
defined in the text. The star that is ejected from the host settles 
into a quasi-stable orbit around the companion. Top panels are 
orbits around the host, bottom panels are around the companion. 
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Figure 8. The change in orbital parameters of a single planet 
around the companion with eo = 0, ao = 10 au, measured at 10 
Myr. These runs were with a 250 au binary and Mi = M2 = 1.0 
M . 



of 10 au, and zero eccentricity. This nearly blank slate of a 
planetary system records the effect of the wandering plan- 
ets on planetary systems near the likely terrestrial and giant 
planet forming region around the companion. 

The changes to the semi-major axis, eccentricity, and 
inclination at 10 Myr of this target planet are plotted in 
Figure [8] The majority of systems are relatively unaltered 
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Figure 9. As in Figure [4] the times spent around each star for 
bouncing systems in a 250 au binary, but with a single planet 
orbiting the companion. Top panels are orbits around the host, 
bottom panels are around the companion. 

by the intruding planets, but of the order 10 per cent of these 
planets experience changes in semi-major axis of around 1 
au, increases in eccentricity to tens of per cent, and inclina- 
tion from a few to a few tens of degrees. 

These changes to the single planet's orbit imply some 
change to the orbit of the intruding planet, possibly altering 
the time spent around each star. In Figure [9] we show these 
times. There is very little change, except for the presence of 
a few systems that spend more time around the companion, 
filling in the space above the diagonal reference line. These 
are not very notable, with two exceptions. Twice the planet 
ejected from the host displaces the planet around the com- 
panion, taking its place and ejecting the companion's orig- 
inal planet from the system. These are the only two cases 
where the planet originally around the companion is ejected. 
One of these runs is shown in Figure 1101 The interactions 
between the incoming planet and the companion's planet 
occur over about 0.5 Myr, culminating in the escape of that 
planet after briefly crossing over to the host system. 

3.6 Non-zero binary eccentricity 

For simplicity, we have dealt exclusively with zero eccentric- 
ity binaries in this first work. The ecc entricity distribution o f 
long-period binaries spans all values; iRaghavan et al.l (| 20101 ) 
find a flat distribution up to e ~ 0.6, with the drop off above 
that value probably due at least in part to observational 
bias. Eccentricity in the binary may expected to change the 
results of this work; Hill stability of a test particle is a func- 
tion of the massive bodies' ec c entricity ijMarchal fc Bozisl 
1 19821 ; iHolman fc Wiegertl Il999l ; ISzenkovits fc Makol l200sK 
and this dependence could alter the timescale or qualita- 
tive nature of the scattered planet's trajectory in the binary. 
While will we defer a more complete study of the effect of 
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Figure 10. Switching followed by planet displacement in a 250 
au, equal-mass binary with a single planet in orbit around the 
companion. The orbital elements are plotted relative to whichever 
star's space the planet is found in at any time, as defined in the 
text. The star that is ejected from the host interacts with a planet 
originally around the companion, before ejecting it and settling 
into a stable orbit around the companion. The oscillations in that 
planet's eccentricity are due to interactions with the host star. 

eccentric binaries on these results to other work, the results 
of one example case hint at the magnitude and character of 
the differences to expect. 

We re-ran the 250 au equal mass case with a bi- 
nary eccentricity of 0.75. The first point to note is that 
the closer approaches of the stars perturbed the planetary 
systems so that only 14 cases remained unscattered af- 
ter the_T^_^Tyj^im^This_destabilizing effe ct is described 
in IHolman fc Wiegertl l| 19991 ) and noted in iMarzari et al.l 
|2005h . Similarly to the latter authors, we find more planets 
ejected (per scattered system) at higher eccentricity; ~ 1.25 
per system compared to ~ 0.95 in the zero eccentricity case. 
The fraction of ejected planets that cross LI is similar, 89 
per cent with e = 0.75. The number of bouncing planets is 
virtually the same at 76 per cent compared to 77 per cent in 
the circular runs. The times spent around the host and com- 
panion are smaller, by about a factor of 2. The introduction 
of moderate eccentricity appears to preserve the qualitative 
and many quantitative aspects of this work, though the de- 
tailed dependence of these results on the binary eccentric- 
ity deserves further attention. We note that as the binary 
periastron decreases to be comparable to the planet's semi- 
major axis the companion may tak e a more direct role in 
the scattering ( Ma rzari et al.l 120051 ) , perhaps leading to a 
qualitative shift in the outcomes. 



4 DISCUSSION AND FUTURE DIRECTIONS 

A natural extension to this work would be to remove var- 
ious restrictions of the CR3BP. For example, if the two 
stars are not bound to one another, but are hyperboli- 
cally scattering off of each other, then one may utilize 
the parabolic or hyperbolic restr icted three-body problem 
l|Luk'vanovll2oTol ; ISari et alJlioioh in order to determine the 
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resulting motion of the planet. However, as touched upon 
in section 13.61 a closer extension to our study would in- 
volve removing the circular orbit restriction but keeping the 
stars bound to each other. If the stellar masses orbit one 
another on an unchanging, eccentric orbit, then the Jaco- 
bian int egral no longer exists but the 5 Lagrangian points 
still do l|Szebehelv fc Giacaglial Hiii ISzebehelvlll967l ). The 
coordinate system traditionally adopted for this situation 
is nonuniformly rotating and pulsating; the zero-velocity 
surfaces change with time. The consequence is a system 
which is diffic ult to describe in simple closed analyti- 
cal forms (e . g. IContopoulosI 1 19671; iDelva fc Dvorak! 1 19791 ; 



IGawlik et alj|2009l; IHou fc Liull201ll) but can still admit es 



cape and capture (e. g. Bailevll 19721; iHuang fc InnanerJI 19831 : 



iLlibre fc Pinolll990l ; lMak6 fc SzenkovitsfeOO^ 



In all of these situations, a star cannot capture a planet 
and retain it over the star's main sequence life without 
some additional dissipative force. In the trad itional pho- 
togra vitational restricted three-body problem (|Radzievskiil 
1950), this force arises from the star's radiation on a mass- 
less dust particle. For a planet in a binary system, however, 
a non- negligible force might arise from additional planets 
(as briefly discussed in Section j3.5l) or from mass loss or 



gain f rom the stars or planet (|Bekovlll988l . Il99ll ; iLuk'vanovl 



2009:; iLetelier fc da Silval 1201 1\ ). This situation 



is particu- 
larly relevant after one or both of the stars have left the 
main sequence, however the bouncing timescales we find here 
(typically less than a few Myr) are shorter than stellar evo- 
lution timescales. For a discussion of planetary transfer in 
post main sequence binary systems, see Kratter & Perets 
(2012, in preparation). 

A more likely dissipative mechanism in the planet 
formation context is interactions with discs, either of 
protoplanets, debris, or gas. Scattering can occur when 
significant amounts of g as are still around the star 



l|Moeckel fc Armitagell2012 '). and if the two stars of the bi- 
nary are at similar stages in their lives a planet scattered 
over to the companion star could encounter a significantly 
dissipative medium. If dissipation is sufficient to strand the 
marauding planet around the companion star, the orbit 
would likely be a wide one, offering a way to produce gi- 
ant planets in the inner regions of a disc and send them off 
to a wider orbit around a binary companion. This sort of sce- 
nario is probably most efficiently explored via a dissipative 
prescription of some sort in the CR3BP rather than direct 
integration of a hydrodynamic disc. Finally, while we have 
focused on planets themselves scattering across to the com- 
panion star in this work, other scattering events are likely 
in planet formation. Scattered planetesimals and exoplan- 
etary analogues to our own outer Solar System structures 
(e.g. the Kuiper belt, or scattered disc objects) may be able 
to participate in jumping between stars in a similar fashion 
to what we have shown here. 

The results of this work are arguably entertaining, but 
are they observable? The short timescales of active bounc- 
ing that we observe would suggest not; the consequences of 
planets bouncing back and forth are likely to be seen in the 
extensions discussed above, particularly a dissipative envi- 
ronment around the companion, which could lead to per- 
manent capture of the scattered planet, or the introduc- 
tion of a planetary system around both stars. The cursory 
examination of the response of a planet around the com- 



panion that we performed here shows that a small but non- 
negligible percentage of compact planetary systems could be 
perturbed by the intrusion of a scattered planet. There are 
many possible consequences of this: changes to the eccen- 
tricity of outer planets can propaga te inwards through the 
system (jZakamska fc Tremaineil2004l ); two-planet systems in 
formally Hill stable configurations could be pushed into un- 
stable architectures; mean motion resonances could be bro- 
ken; and richer systems than the single planet that we simu- 
lated could be induced into scattering themselves, resulting 
in wholesale transfer of planets between the two stars. These 
are the topics of ongoing investigations. 



5 CONCLUSIONS 

Via direct integration of a planetary system around one 
member of a circular binary system, we have examined the 
idea that planets can be sent bouncing between the stars 
during planet-planet scattering. The internal dynamics of 
the planetary system can increase the energy of one or more 
planets, enabling it to cross over to the companion star from 
its host. We find that this is in fact a dominant outcome of 
scattering in a binary: about 70 to 85 percent of planets 
ejected from binaries in the separation range 250 to 1000 au 
will spend time around both stars before leaving the system 
entirely, and about 45 to 75 percent will bounce repeatedly 
back and forth. The amount of time spent bouncing is as- 
tronomically short: a few to tens of binary periods, although 
with some systems surviving hundreds to thousands of peri- 
ods. We propose that the impact of this bouncing will most 
likely be seen in the response of the companion star's own 
planetary system, or with the inclusion of un-modeled dis- 
sipative forces, such as discs, around the companion. 
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